Zero bias anomaly in a two dimensional granular insulator 
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We compare tunneling density of states (TDOS) into two ultrathin Ag films, one uniform and one 
granular, for different degrees of disorder. The uniform film shows a crossover from Altshuler- Aronov 
(AA) zero bias anomaly to Efros Shklovskii (ES) like Coulomb gap as the disorder is increased. 
The granular film, on the other hand, exhibits AA behavior even deeply in the insulating regime. 
We analyze the data and find that granularity introduces a new regime for the TDOS. While the 
conductivity is dominated by hopping between clusters of grains and is thus insulating, the TDOS 
probes the properties of an individual cluster which is "metallic". 
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The metal-insulator transition (MIT) in disordered 
electronic systems renders their transport properties one 
of the most exciting topics in condensed matter physics. 
The effects of Coulomb interaction in the vicinity of the 
MIT are particularly strong, and despite major theo- 
retical and experimental efforts are not yet fully under- 
stood [2j. Much of the work is focused on 2D metallic 
films. While infinite 2D systems are insulators [lj, sam- 
ples shorter than the localization length £ exhibits metal- 
lic behavior. This enables to utilize thin dirty metals to 
study both metallic and insulating regimes. 

For weakly disordered metallic films with relatively 
high conductivity the interplay between electron-electron 
interaction and disorder has a pronounced influence on 
the electronic properties, as shown by Altshuler and 
Aronov (AA) [3|. Both conductivity and the tunneling 
density of states (TDOS) are significantly modified (rel- 
atively to the non-interacting values) and acquire tem- 
perature dependent corrections. The theory predicts 
the depletion of the TDOS around the Fermi level, a 
phenomenon that became known as zero bias anomaly 
(ZBA), which was indeed observed experimentally in dis- 
ordered metals |J| . 

For strongly disordered systems, electronic states are 
localized, and transport at low temperatures is achieved 
via variable range hopping (VRH). For non- interacting 
electrons there is no gap in the TDOS at the Fermi level 
and the electric conductivity is given by Mott's law [5|. 
Electron interactions open a Coulomb gap at the Fermi 
energy, affecting the TDOS and the electric conductivity 
as shown by Efros and Shklovskii (ES) [6(. The Coulomb 
gap in the 2D density of states is given by 
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and the conductivity as a function of temperature is 
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FIG. 1: Schematic description of the system under consider- 
ation: Metallic grains of the size L form clusters with typical 
size £. The dashed lines denote the weak links, connecting 
different clusters. 



where k is a dielectric constant and To 



7«£. 




Though the ZBA and the Coulomb gap both result 
from the long range Coulomb interaction and lead to a 
depletion of the TDOS at the Fermi energy, a unified 
theory of both effects is still absent [8[ . 

Most theoretical and experimental works deal with ho- 
mogeneously disordered systems in which the degree of 
disorder controls the vicinity to the MIT. A different sys- 
tem that may provide a natural and convenient arena for 
experimental and theoretical study of strongly interact- 
ing systems is a granular metal |9| i.e. a disordered array 
of metallic islands with typical size L, embedded in an 
insulating matrix (schematically shown in Fig. [1}. In 
these systems the conductivity of the grains, <7 is high, 
so that conductivity of a sample is determined by the 
inter-grain coupling g s , which corresponds to the Drude 
conductivity a = g g e 2 /h, and consequently to the diffu- 
sion coefficient D g = g s /vo- Clearly, the intergrain D g 
and g g are very different from the metallic go and Dq 
within the grain. 

Due to the randomness in the position of the grains 
the coupling strengths g g fluctuates. Deep in the insu- 
lating regime the transport between the grains is gov- 
erned by interaction-controlled VRH resulting in the ES 
formula, Eq.©. Typically £ > L so one may visual- 
ize the sample as a network of weakly connected clusters 
(with cluster size £). Each cluster contains several grains, 



and is connected to a neighboring cluster by weak links 
with dimensionlcss conductivity g c . The measured con- 
ductivity a(u),q), is controlled by either g c <C g g <C go, 
depending on the wave vector q. In the dc limit the 
conductivity of the sample is determined by g c . The 
TDOS, on the other hand, probes the local properties 
of the system, and as a function of energy explores sev- 
eral regimes, as is illustrated in Fig. 2. At energies larger 
than 62 = e 2 <7 g /27rC (C ~ 2ttkL is an effective capaci- 
tance of a grain; tac = ^/^2 is the corresponding charg- 
ing time of the grain) the tunneling probes the properties 
of the system on spatial scales smaller than the size of a 
grain. Depending on grain's geometry, the TDOS in this 
case is the same as in the bulk (2D or 3D) metal and is 
described by perturbativc AA theory. For energies below 
61 = -D c /£ 2 the TDOS explores the universal low energy 
physics. In this limit the granular structure of the ma- 
terial is irrelevant and one recovers the Coulomb gap of 
Eq.© in the TDOS. 

For the intermediate energy range ei <C e <C £2 one 
probes spatial distances shorter than the cluster size £, 
but larger than the size of a single grain. This is a new 
regime, that is unique to the granular material. While 
the dc conductivity at this energy regime is governed by 
tunneling between different clusters, the TDOS is deter- 
mined by the properties of a single cluster. The tunneling 
in this regime occurs from the reference electrode into a 
grain, which is strongly coupled to other grains in the 
same cluster. The later can be viewed as the environ- 
ment. 

This kind of problem has been extensively studied in 
the context of the environmental ZBA in nanostructures 
[7(. It was shown that as far as TDOS is concerned, the 
original problem can be mapped onto an effective elec- 
tric circuit. In this description, the microscopic details 
of the problem arc encoded through the impedance of 
the circuit. In accordance with fluctuation-dissipation 
theorem this impedance determines the spectrum of the 
fluctuations of the electromagnetic field. Accounting for 
these fluctuation one calculates the ZBA that develops on 
time scales ranging between tuc (the time it takes the 
system to screen an added electron) and t\ = %/ei (the 
time at which the charged cloud reaches the boundaries 
of the conducting cluster). From this consideration, it is 
clear that in this intermediate energy range the problem 
is completely equivalent to the one solved in the context 
of transport through an open quantum dot [1CH14J . The 
tunneling conductance between the dot (analogous to the 
grain) and the lead (analogous to the cluster) yields a 
ZBA as a function of temperature and voltage (see Ref. 
12j for the details). To compare the theoretical predic- 
tions with the experimental observation we analyze the 
exact results (written in terms of the digamma function) 
in various asymptotic limits. Inside the narrow energy 
interval e\ < e < gA/2ir (of the order of IK, see esti- 
mates below) the TDOS only weakly depends on energy 
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FIG. 2: A schematic description of the dependence of the 
TDOS on energy for a granular system. In the high energy 
limit (e > £2) one observes the TDOS of a disordered metal, 
accounted for by the Altshuler-Aronov theory. In the low 
energy limit (e < 61) the system is in the VRH regime and 
exhibits a Coulomb gap given by Eq.(fl]). For intermediate 
energies (ei < e < £2) the electric conductivity is given by 
Eq. ([2]) while the TDOS is determined by tunneling into a 
cluster, Eq.®. Insert: TDOS, given by Eq.JTJl, for e < ei. 

and may be approximated by a constant; here A is the 
energy level spacing between non-interacting single par- 
ticle levels in the dot. 

For the experimentally relevant energy range gA/2ir -C 
eCf2 one obtains logarithmic ZBA 



v(e)/v = H ln(etflc) 
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The applicability of this result for tunneling into a sin- 
gle open quantum dot, coupled to the environment, was 
recently confirmed experimentally in Ref. [la ]. Note, 
that Eq.Q can be also viewed as the zero dimensional 
version of the AA theory for short range interaction. 
The presence of an external metallic tunneling electrode 
makes this model rather realistic. 

In this Letter we compare the TDOS of two thin dis- 
ordered films, uniform Ag and granular Ag, as a function 
of g and experimentally detect a regime in which Eq.Q 
determines the TDOS of the granular film. In this unique 
regime the TDOS depends logarithmically on energy de- 
spite the fact that g extracted from conductivity is much 
smaller than unity. 

The disordered 2D samples were prepared by quench 
condensation i.e. sequential evaporation of ultrathin lay- 
ers of Ag on a cryo-cooled substrate within the measure- 
ment apparatus. This technique enables one to study the 
transition from insulating to metallic behavior on a sin- 
gle sample without thermally cycling it or exposing it to 
atmosphere. Depending on the choice of substrate, one 
can fabricate either granular or uniform disordered films. 
If the samples are quench condensed on a bare SiO sub- 
strate, the resultant film is discontinuous containing sep- 
arated Ag islands about lOnm in diameter [lfj. Adding 



material results in a decrease of the average distance be- 
tween the islands leading to an exponential drop of the 
resistance with thickness. On the other hand, when the 
substrate is pre-coated by a thin layer of amorphous Ge or 
Sb the deposited film is electrically continuous at a thick- 
ness of 1-2 monolayers of material (~5A) [l7|, [l8| after 
which the Ag film obeys the usual Ra ocl/d dependence, 
where d is the sample thickness. These films are referred 
to as " uniform films" . We note that one cannot entirely 
rule out some granularity in these films as well, however, 
the fact that finite conductivity is achieved for a mono- 
layer of evaporated material followed by ohmic behavior 
implies that the samples are very homogeneous. 

For measuring TDOS we fabricated tunnel junctions 
using the following scheme: Prior to the cool down, four 
gold leads were thermally evaporated on an insulating 
Si/Si02 substrate. Then a 25nm thick Al strip with di- 
mensions of 1.5mm by 5mm was e-beam evaporated so 
that it connected two of the gold leads. The Al strip 
was then exposed to atmosphere for 2 hours thus form- 
ing an oxide layer of AI2O3 with a thickness of 3 — 5nm 
that served as a tunneling barrier. The substrate was 
then connected to the Hc3 pot of a He3 fridge which was 
pumped to ultra high vacuum to enable film evapora- 
tion. Transport measurements were performed on the 
Al layer to confirm the presence of superconductivity 
with Tc ~ 1.2K. Then a 1.5mm by 2mm Ag strip 
was quench condensed cither directly on the substrate 
for a granular film or on a pre-evaporated 2nm thick 
layer of Sb for a uniform film, thus connecting the other 
two gold leads as shown in the inset of Fig. [3^- Deposi- 
tion rate and film thickness were monitored, in situ, by a 
quartz-crystal. The evaporation was terminated at any 
desired resistance. After each step of growth an R(T) 
measurement of the Ag film determined the degree of 
disorder. Then, the dl/dV versus V across the junctions 
were measured using standard lock-in techniques while 
making sure that Rbamer >> Rsampie so that the Ag 
film could be regarded as an equipotential electrode even 
in the insulating state. All measurements were performed 
at a temperature below Tc of the Al. We first made sure 
that a BSC gap could be detected and then applied a 
5KG magnetic field to suppress superconductivity and 
measure the TDOS of the disordered film. 

The differential conductance, dl/dV, representing the 
TDOS, as a function of voltage for both films is shown in 
Fig. 03 for different values of g. The uniform films shows a 
transition from a dl/dV oc ln(V) dependence for g = 13 
(weak disorder) to dl/dV oc V for g = 0.35 (strong dis- 
order). This could be interpreted as a crossover from 
AA ZBA to an ES Coulomb-gap-like behavior. Indeed 
the AA theory is applicable only for j > 1, while for 
g < 1 one may expect strong localization physics to dom- 
inate. A similar crossover was reported by Butko et-al 
on thin films of Be [19j. It should be noted, however, 
that probing the Coulomb gap by a tunnel junction is 




FIG. 3: dl/dV versus V of uniform film having g = 13, 1, 0.35 
from top to bottom respectively (a) and a granular film having 
g = 32, 8.6, 3.25, 0.2 from top to bottom respectively (b). The 
solid lines are fits to Eq. [3] The inset in (a) shows a schematic 
description of the junction with a granular Ag film as the 
top electrode. The inset of (b) depicts the resistance versus 
temperature of the most insulating stage of the granular film. 
The solid red line is a fit to the ES law of Eg. ^ (color 
online). 

problematic. On the one hand the insulating layer must 
be thin to allow for tunneling, but on the other hand the 
close electrode screens the Coulomb interaction over dis- 
tances larger than the thickness of the tunneling barrier. 
This severely modifies the Coulomb gap which originates 
from long range interaction. It seems therefore that while 
tunneling experiments are a good tool to probe the im- 
portance of Coulomb interactions in Anderson insulators, 
they do not yield accurately the shape of the Coulomb 
gap. 

The situation in the granular film is different. Fig. |3Jd 
shows that the dl/dV curves follow a ln(V) throughout 
the entire crossover from weak to strong disorder. This 
means that even deep into the insulating regime, where g 
extracted from conductivity is as low as 0.2 the TDOS of 
the granular film exhibits AA physics. The inset in Fig. 
[3Jd shows that the conductivity in this film follows the 
ES law of Eq.©, so naively one would expect the TDOS 
to exhibit a curve more similar to the ES Coulomb gap, 
Eq. (TTJ) . Clearly, this is not the case for the granular film. 
The solid lines in Fig. [3] arc fits to Eq.Q. It is seen 
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FIG. 4: The dimensionless conductance, g, and the charging 
time tnc as a function of Ra in units of h/e 2 . Symbols are 
results obtained from the fits of the data in Fig. [3] (a) to Eq. 
[3] The dashed line in the left panel is the expected value from 
g — 1/Rn ■ The dashed line in the right panel is the calculated 
value assuming a lOnm grain. 

that all the granular sample curves fit very well to this 
expression. In contrast, it was impossible to obtain a 
reasonable fit for the curve of the insulating uniform film 
(g = 0.35). 

These fits enable one to extract the charging time tnc, 
and the conductance g as a function of the disorder. 
Fig. [4] depicts these two fitting parameters as a func- 
tion of R a ■ The extracted tfic increases monotonically 
with R a as expected and is very close to that calculated 
value from Eq. Q based on the capacitance of a lOnm 
dot. This reflects the fact that, with increasing distance 
between grains, the electron spends more time in the 
cluster. On the other hand g obtained from this fit is 
much larger than that extracted from the conductivity 
g = -s-\- Clearly, conductivity and TDOS probe a dif- 
ferent g. While conductivity is sensitive to the weak cou- 
pling between clusters, g c , the TDOS apparently probes 
g g which corresponds to the electrical coupling between 
the grains within each cluster. 

Returning to the sample morphology described in Fig. 
[I] we note that a typical size of a grain is L ~ lOnm 
|20( while £ extracted from the R(T) curve of the inset 
of Fig. [SJd utilizing Eq.© is found to be ~ lOOnm. This 
means that £ 3> L and the analysis presented above is 
appropriate. Using g g ~ 10 for metallic conductivity we 
get for our samples: g/S.jl-K ~ 0.25meU and ge 2 /2ttkC ~ 
30 — lOOmeV. Hence Eq.Q should hold for the gate 
voltage regime: 0.3mU < V < 30mU . The experimental 
region of the gate voltage were chosen within this range, 
thus the TDOS is indeed expected to probe the physics 
of a cluster of grains and is should not be sensitive to the 
inter cluster weak connections. 

To summarize, we have studied the electric conductiv- 
ity and the TDOS in homogeneous and granular films 
for various strengths of disorder. The homogeneous sam- 
ple is characterized by a single energy scale, t\ ~ D/^ 2 , 
above which the TDOS exhibits a metallic ZBA and be- 
low which there is a Coulomb gap. This energy becomes 
rather large as the sample is driven deep into the insu- 



lating regime. Hence, a logarithmic dependence of the 
TDOS on energy will not be observed experimentally 
in strongly disordered homogeneous films. A granular 
structure gives rise to a new energy interval, in which 
ZBA is similar to the environmental Coulomb blockade 
in nano-structures. This results in logarithmic depen- 
dence of TDOS on e down to relatively low energies and 
explains the observation of AA type ZBA in a highly 
disordered system. 
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